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Abstract: Boost-invariant dynamics of a strongly-coupled conformal plasma is studied 
in the regime of early proper-time using the AdS/CFT correspondence. It is shown, in 
contrast with the late-time expansion, that a scaling solution does not exist. The boundary 
dynamics in this regime depends on initial conditions encoded in the bulk behavior of 
a Fefferman-Graham metric coefficient at initial proper-time. The relation between the 
early-time expansion of the energy density and initial conditions in the bulk of AdS is 
provided. As a general result it is proven that a singularity of some metric coefficient 
in Fefferman-Graham frame exists at all times. Requiring that this singularity at r = 
is a mere coordinate singularity without the curvature blow-up gives constraints on the 
possible boundary dynamics. Using a simple Fade resummation for solutions satisfying 
the regularity constraint, the features of a transition to local equilibrium, and thus to the 
hydrodynamical late-time regime, have been observed. The impact of this study on the 
problem of thermalization is discussed. 



Keywords: Gauge-gravity correspondence. Heavy lon^ . 



'guillaume .beuf@cea.fr, michal .hellerOuj . edu.pl , uf rjanikOth. if .uj . edu.pl , 
robi . peschanskiOcea . f r 



Contents 



|l]. Introduction ^ 

|2]. Boost-invariant dynamics and holography ^ 

pTl| Energy-momentum tensor of boost-invariant flow ^ 

T2| Holographic reconstruction of space-time S 



|3|. Late time dynamics ^ 

p] Motivation | 

3.2 The asymptotic solution ^ 

3.3 Large proper-time expansion and hydrodynamics ^ 

^ Early time dynamics 10 

[4.1| No scaling at early time 10 



4.2 The Feffer man- Graham metric is singular at all times 12 
|4.3| Constraints at early time 14 
Early time expansion of the energy density 16 



|5]. Transition to the hydrodynamic regime 17 
|5.1| Resummation scheme for the energy density 17 
Qualitative features of the approach to local equilibrium ff9 



IgL Summary 20 



1. Introduction 

Recently, a new theoretical approach to heavy-ion reactions and the formation of quark- 
gluon plasma (QGP) in RHIC experiments using the AdS/CFT correspondence arose as a 
fruitful application of string theory to the real world. The AdS/CFT correspondence allows 
to describe many properties of a class of gauge theories at strong coupling in terms of a 
dual gravitational description. What makes it a unique tool is firstly that it works well in 
Minkowski signature allowing for a description of inherently time-dependent phenomena, in 
contrast to lattice QCD methods which are directly tied to Euclidean signature. Secondly, 
explicit analytical computations are often possible and one can get new insight into strongly 
coupled dynamics from the dual geometric perspective. 

On the experimental side, the features of the distribution of particles observed in the 
final stage of the reaction point to the existence of a hydrodynamical regime of the QGP 
expansion |2|, in particular the observed elliptic flow denoting sizeable collective effects 
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hydrodynamics 



Figure 1: Description of QGP formation in heavy ion collisions. The kinematic landscape is 
defined hy t — (x^Y ~ (x^)'^ ; y = ^ log fn^fj ; x±—{x^,x^} , where the coordinates along the 
light-cone are x° ± cc^, the transverse ones are {x^, x^} and r is the proper time, y the "space- time 
rapidity" . 



Q. The hydrodynamic regime has to last long enough and start soon enough after the 
collision in order to explain the observed collective effects. Moreover, the smallness of the 
viscosity which can be extracted from hydro dynamical simulations describing the data leads 
to an almost-perfect fluid behaviour of the QGP, and thus to a short mean-free path inside 
the fluid. Putting together these experimental inputs, and in order to go beyond a mererly 
phenomenological description, it appears to be theoretically necessary to investigate as 
much as possible the properties of a strongly-coupled Quantum-Chromodynamic plasma. 

In the absence of nonperturbative methods applicable to real-time dynamics of strongly 
coupled Quantum Chromodynamic (QCD) plasma, one is led to consider similar problems 
from the point-of-view of the AdS/CFT correspondence, that is looking for the characteris- 
tics of plasma in a gauge theory for which the AdS / CFT correspondence takes its simplest 
form - the Af = 4 supersymmetric Yang-Mills theory [Q| which possesses a known and 
tractable gravity dual. 

Although the AA = 4 gauge theory is supersymmetric and conformal and thus quite 
different from QCD at zero temperature, both supersymmetry and scale-invariance are 
broken explicitly at finite temperature and we may expect qualitative similarities with 
QCD plasma for a range of temperatures above the QCD deconfinement phase transition^. 

Indeed, the gauge/gravity dual calculation ^ showing, in a static setting, that the 
viscosity over entropy ratio r]/s is very small (equal to l/47r) and even suggesting a universal 
lower bound, is in qualitative agreement with hydrodynamic simulations of QCD plasma 
and was a poweful incentive to explore further the AdS / CFT duality approach. 

In order to go beyond static calculations, one has to adapt the dual AdS/CFT approach 
to the relativistic kinematic framework of heavy-ion reactions, where two ultra-relativistic 
heavy nuclei collide and form an expanding medium, see Fig.|^. It is convenient, initially, 

^ There exist more refined versions of tlie AdS/CFT correspondence wfiich may fiave more features in 
common witli QCD, however the gravity backgrounds are much more complicated and we will not consider 
them here. 
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to adopt the simplest setting in which the relevant physics is still present. 

For this sake, one is guided by the seminal work of Bjorken who made a simplifying 
assumption of boost-invariant dynamics for his hydrodynamic description of high-energy 
heavy-ion collisions. Hence, the different stages of the plasma formation, expansion and 
hadronization, see Fig.|l|, depend only on the proper-time r and not on rapidity^ y. This 
assumption is indeed phenomenologically valid for the central rapidity region of the collision 
at high energy. 

For the application of the AdS/CFT correspondence to this setting, and in particular 
for finding the gravity dual |^ , the above symmetry assumptions allow for a drastic simpli- 
fication of an otherwise yet probably intractable problem, since the dual dynamics depends 
only on two variables, r, proper-time and z, the bulk variable in the dual asymptotically 
AdS spacetime. 

The two main physical questions^ are to understand the behaviour of quark gluon 
plasma at late and early times. At late times, we would like to derive the appearance of 
hydrodynamic behaviour with all its characteristics, asymptotically perfect fluid flow, low 
viscosity, transport coefficients, etc. At early times, we would like to understand the rapid 
thermalization starting from the ultra-relativistic initial conditions. 

In fact the study of the former problem has led to quite successful developments, 
showing that the hydrodynamic boost-invariant Bjorken flow is a consequence of late-time 
dynamics of gravity in the dual side of the AdS/CFT correspondence. Subsequently the 
way in which hydrodynamics arises from the AdS/CFT has been understood even without 
using any symmetry arguments [^j as a gradient expansion of Einstein's equations. However 
since our main interest in this paper is the radically different early time regime where such 
expansions are inapplicable, we will use the full system of Einstein's equations with the 
boost-invariant symmetry assumptions. Before that, we will first recall the analysis for late 
times, since the point of departure will be analogous for the early time problem, but with 
different developments and results. 

The second problem, which forms the novel part of our investigations in the present 
study, is to go beyond hydrodynamics and investigate early-time dynamics. The only 
assumption we will make, apart from boost-invariance, is the existence of a strongly coupled 
pre- hydrodynamic stage^ allowing for the use of the AdS/CFT correspondence in the same 
framework as for late times, and thus to look for a connection of the two regimes in a 
unified scheme. In particular, one is curious to understand the transition towards the 
hydrodynamic regime and thus to get some information about the puzzling problem of the 
apparently short thermalization time. 

Summarizing the contents of the paper, the next section recalls the basic AdS/CFT 
approach to boost-invariant dynamics. The main results of the late-time solution are 

^For the Bjorken flow, spatial and momentum rapidity are equal. 

^We do not address here another important problem which is hadronization, since it requires to go 
beyond our current understanding of the presently known Gauge/Gravity dualities. 

^Physically, it may be that the initial stage of the collision has some weak-coupling properties, for 
instance in the Color Glass Condensate (CGC) setting but we will not consider here this possibility (for 
a recent attempt to describe CGC initial conditions within the AdS/CFT see ]ic|]). 
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given in section 3, with emphasis on those aspects which will be useful for the early-time 
regime. Section 4 is devoted to the theoretical analysis of the early-time boost-invariant 
dual dynamics, while section 5 deals with numerical solutions with the goal of studying the 
transition to hydrodynamics. Conclusions and outlook are given in the final section. 



2. Boost-invariant dynamics and holography 

2.1 Energy-momentum tensor of boost-invariant flow 

Currently the best understanding of non-perturbative non-equillibrium 4-dimensional quan- 
tum field theories is provided by string theory methods. Holography translates gauge theory 
problems in four dimensions at large number of colors and strong ('t Hooft) coupling into 
higher-dimensional classical gravity. Einstein's equations involve then, in the simplest set- 
ting, functions of both gauge theory dimensions and of one more variable z parameterizing 
the radial direction in AdS. In a general setup, finding solutions of Einstein's equations 
must involve numerical analysis and is rather difficult. Moreover, singular components of 
the metric may and indeed do appear ||^, 12] which are not easy to implement numer- 



ically, unless their structure is known beforehand. It is thus important to focus on the 
simplest possible, yet phenomenologically interesting, gauge theory dynamical processes. 
As explained in the introduction, we will focus on the boost-invariant plasma expansion. 
This setup serves as a first approximation of dynamical plasma created at RHIC, as ex- 
cited nuclear matter undergoes one-dimensional expansion and is both translationally and 
rotationally invariant in the plane perpendicular to the expansion axis (physical quantities 
do not depend on The crucial simplifying assumption is boost-invariance along 

the expansion axis, which can be most easily imposed when passing to proper-time (r) 
and (spatial) rapidity (y) coordinates related to the usual lab-frame time (x^^ and position 
along the expansion axis (x?) by 



X = rcoshy 

x^ = rsinhy . (2-1) 

Under a boost along the x^ direction proper-time does not change, whereas rapidity is 
shifted by a constant. Thus in this coordinate frame boost-invariance means that physical 
observables do not depend on rapidity. 

A quantity of special interest is the expectation value of the energy-momentum tensor 
of the plasma. In the (r, y, x^, x^) frame with all mentioned symmetries^ imposed, the most 
general boost-invariant energy-momentum tensor takes the form 

V = diag {e (r) , t\\ (t) , (r) , (r)} . (2.2) 

Imposing energy-momentum conservation V^T'^'^ = (in boost-invariant case there is only 

one independent equation) and tracelessness T^* = leads to the energy momentum tensor 

fully specified in terms of a single function of one variable - the energy density e (r) , namely 

^Together with parity invariance in rapidity, which ensures that there are no nondiagonal Try terms in 
the energy-momentum tensor. 
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T,, = diag (r) , -r\e (r) + re' (r)), e (r) + ire' (r) , e (r) + ^re' (r)} . (2, 



3) 



The function e (r) may be arbitrary at this level and it is the gravitational dual which 
will single out the allowed ones. Note that T^y is submitted to an important positivity 
constraint, namely T^^t^^t^ > 0, for any time-like vector implying 0] 

< < 0. (2.4) 

r 

This requirement may perhaps be lifted in case of some highly quantum states |jl^ ap- 
pearing at transient times, but we expect that it should hold for any more macroscopic 
configuration. We observe a temporary violation of this constraint for a specific initial 
condition (see section 5). 

2.2 Holographic reconstruction of space-time 

The AdS/CFT correspondence maps the expectation value of the energy-momentum tensor 
operator into the behavior of the 5-dimensional asymptotically AdS metric Gab- This 
metric is a solution of Einstein's equations with negative cosmological constant A = — 6L~^ 
(where L is the AdS radius^), namely 

Rab -\r-Gab-Q Gab = 0. (2.5) 

The general asymptotically AdS metric written in Fefferman - Graham coordinates ||I 
takes the form 

d«^ = GABdy^dy^ = gM^dxMx- + dz^ 

where z is the AdS radial variable and g^^^j is a four-dimensional metric on constant z slices, 
which depends both on the field theory directions and z. The hypersurface at z = is a 
conformal boundary of AdS and is identified with the metric in which the CFT lives. In the 
considered field theory setup this is the 4-dimensional Minkowski metric ry^i^. Einstein's 
equations obtained from (|2.5| )) (see further Eq. ( [4.10 ) for the full set of Einstein's equations 



expressed in Fefferman- Graham coordinates) can be solved pertubatively in the vicinity of 
the conformal theory. The four-dimensional part g^^ of the AdS metric'' 

9,. = af} + 5g + 9^} + 9f} + +9^} z' + ... (2.7) 

can be obtained order by order in z'^ from the equations of motion ( [2.51 ). The leading piece 
in ( ^.7] ) is the Minkowski metric r/^jy, which in the proper (r, y, x±) variables takes the form 



®In the rest of the paper the AdS radius L will be set to 1, but can be trivially restored on dimensional 
grounds. Greek and latin indices denote respectively field theory and bulk AdS directions. 

^AU the formulas regarding holographic reconstruction of the bulk space time are provided for the case 
when the boundary is the Minkowski space. 



- 5 - 



gl^Jdx^'dx'' = 7]^^dx^'dx'' = -dr^ + r^dy^ + dx^ . (2.8) 

The quadratic term gj^J vanishes identically, and the first non-trivial term g^i^J turns out 
to be proportional to the expectation value of the energy-momentum tensor [^] 

#J = ^ {T,u) . (2.9) 

Higher order contributions [gj^J^^ for n > 2) can be obtained from the Einstein's equations 
and are fully expressed in terms of {Tf^u) and its derivatives. At this level there are no 
constraints on the form of the energy-momentum tensor - any conserved and traceless one 



will do (apart from the positivity conditions (2.4)). 



However it is to be expected [0] that generic energy-momentum tensor would lead to 
singularities in the bulk of AdS, which are beyond the scope of a perturbative solution in 
z^. Therefore in order to have a control over the geometry in the region where potential 
singularities may develop, one needs to solve the Einstein's equations deeper in the bulk. 
Since the energy-momentum tensor for the boost-invariant flow takes a particularly simple 
diagonal form ( |2.2D , in this case we are led to adopt the 5-dimensional Fefferman-Graham 
metric ansatz satisfying the symmetries of boundary dynamics 

_e«(^.^)dT2 + T2e''(^'^)dy2 + e^^^'^Ux^: + dz^ , , 

ds" = K . (2.10) 

The full set of Einstein's equations corresponding to the parametrized metric (|2.10| ) are 
given in ( 4.10| ). Solving the Einstein's equations for any value of r and z is very difficult and 



must probably involve numerical procedures. However it is possible to gain nevertheless 
some analytical insight. There are two regimes where some simplifications may occur - 
these are for small or large values of r. In both cases there is a small parameter (either r 
or its inverse), which can be used to construct the geometry in perturbative manner in r, 
but exactly in z. So we shall focus on late and early time dynamics. 



3. Late time dynamics 
3.1 Motivation 

At sufficiently late times, the boost-invariant plasma is expected to be locally equilibrated 
and thus well-described by hydrodynamics. However, it is a priori not obvious how to 
prove it on the gravity side of the AdS/CFT. The solution was first shown in and then 
further developed in [|16|, |ll|, |12|, |2^, |2l| (see [||] for a review). The following 



subsections review how to obtain the late-time boost-invariant dynamics perturbatively, 
since the methods developed for this approach will be crucial in the rest of the paper when 
considering the open problem of early time dynamics and the corresponding transition to 
local equilibrium. 
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3.2 The asymptotic solution 

The goal of the whole approach both in the early and late times domain is to solve the 
Einstein's equations ( [4.10I ) perturbatively in proper time r, but (at each order in r) exactly 
in the radial AdS variable z. This can be achieved by resumming the power series in for 



the bulk metric (2.7). For arbitrary proper-time dependence of the energy density it is yet 
out of reach, but focusing on the large proper time regime and assuming that the energy 
density of the boost-invariant plasma exhibits a power-like tail 

6(r)~4+- (3-1) 

for some positive power^ of s this can be basically done. The crucial observation of Ref. Q 
is the existence of a scaling variable v = zt~^/'^, which captures the large proper-time 
behavior on the gravity side (see section ([4.1| ) for more details). Taking the r— >oo limit 
while keeping the scaling variable fixed reduces the Einstein's equations ( [4.10 ) for the 



functions a(T, z), b{T,z) and c{t,z) to ordinary nonlinear differential equations for their 
scaling forms gq (v), bo (v) and cq (v), which can be solved exactly leading to 



oo (v) = A{v) — 2m (v) , 

bo{v) = A{v) + {2s-2)m{v), 

co{v) = A{v) + {2-s)m{v), (3.2) 



where 



Aiv) = i{log(l + A(s) v^)+log{l-A{s) v^)}, 

m (v) = {log (1 + A {s) v^) - log (1 - A (s) v"") } (3.3) 



with 



A„) = ,/5£i^, ,3,4) 



It turns out however, that the geometry is generically singular for all the values of s, apart 
from s = 4/3. This can be checked by evaluating the square of the Riemann tensor (i.e. 
Kretschmann scalar) in the scaling limit and requiring its regularity (the authors of Ref. p3|| 
made also the observation that s = 4/3 is the only value for which the functions (|3.2| ) are 
single- valued). 

For s = 4/3 the energy-momentum tensor takes asymptotically the form 



*The power s is constrained to take values in the range < s < 4 assuming positivity of the energy 



density (2.4) in any time-like reference frame. There are three particular phenomenologically interesting 
values of s within this range: s = 4/3, s = 1 or s = 0. The first one leads to the perfect fluid case where 
Pj_ (r) = p|| (r) and the second one to the free streaming scenario with (r) = 0. The third interesting 
regime discussed in the literature ji^, is the limiting case s = with p± (r) — — (r) which will be 
adressed in section 5. 
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r^, = diag|-|j + ...,r2i-^ + ...,i-|j + ..., + (3.5) 

where ^ denotes some numerical constant. From this expression it is clear that (r) = 
P± (r) (asymptotically), the system is in local equilibrium and is described by the perfect 
fluid hydrodynamics. The gravity dual in this scaling limit is given then by 



(3.6) 

and looks like a boosted and dilated black brane 

(l-z'^X'^S^ 1 1 

= - o,, , .:.. u^ dx^ dx^ + ^ (1 + ^^^^) (^M^' + ^M^^) dx'^dx'^ + ^dz2 (3.7) 
z"^ (1 + z**A'*) / 

with boost parameter = \ ■ [drY" and dilatation A ~ t~^/^ . Recently this connection 
has been exploited and formulated as a fluid/gravity duality [^] . Note that ( p.6| ) contains a 
(tt) coefficient which goes to zero at z ~ t^^^ , giving a hint for a black brane horizon |24]. 



It is also clear that the metric (|3.6| ) is not an exact solution of Einstein's equations - 
there are subleading effects not captured by the scaling variable limit. Those with power- 
like scaling correspond to dissipative corrections in hydrodynamics. 

3.3 Large proper-time expansion and hydrodynamics 

When evaluating 7i? = TZabcdT^^^'^^ in the scaling limit {t ^ oo keeping v = z ■ t^^/^ 
fixed) on the asymptotic solution, the following pattern is encountered 

7^2 = (nonsingular)7?.^Q-) (v) + ■ (singular)7^^4/3-, {v) + . . . (3.8) 

where (singular)7?.^^y3-j (v) has a fourth-order pole at f = 3^^'^ and the dots ". . ." denote 
terms suppressed by inverse powers of r higher than 4/3. Any background with singularities 
not covered by the event horizon would correspond to unphysical configurations on the 
gauge theory side (see e.g. |25, ^). Thus the condition of nonsingularity is the crucial 



requirement, which singles out the correct behavior of the geometry and thus of the energy- 
momentum tensor of the boundary plasma. It is obvious that the singularity at r~^/^ in 
( |3.8| ) can be cancelled only if there is an additional contribution to the metric warp factors 
which arises from a correction to the perfect fluid asymptotics 

e(r) = 4(l + ^ + ---) (3.9) 

Assuming now the expansion 

a (r, z) =ao[z- r-^/^) + ^ a, (^z ■ r-^/^) + -1^ [z ■ r-^/^) + . . . (3.10) 
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with analogous equations for the other warp factors b and c and solving Einstein's equations 
at orders t~'^ and t~'^^^ for r ^ 2/3 one ends up with a singular Kretschmann scalar TZ^ 
of the following form 

TZ^ = (nonsingular)7?.^Qj (?;) H — - ■ (nonsingular)?^^^,-^ (v) + (3-11) 
(singular)7^(4/3^ (v) + ■ (singular)7^(2r•) (^) 



However for r = 2/3 the singular contribution at r"^** cancels the singularities at t~^/^ (see 
12] for details) and gives a regular Kretschmann scalar up to the terms of order r~^/^ in 



the scaling limit. Physically, on the gauge theory side, this exactly corresponds to the first 
viscous corrections to the perfect fluid expansion. This argument repeated at any higher 
order fixes the large proper-time expansion of the bulk metric to take the form 



.-1/3 

1 . / 1 



+ ... 



c(r..) = CO + J^e, (..-V3) + J_,, („-./3) + , .. (3.12) 

The functions ai{v), bi{v) and Ci{v) can be obtained order by order solving Einstein's 
equations in the scaling variable. At each order i > there is a single integration constant 
fixed neither by the equations of motion nor the asymptotic behavior of the metric, which 
is a combination of various transport coefficients of the plasma. Each constant^ in the 
expansion of the energy density or equivalently the temperature (e ~ T^) 

^' t'/^X Ar2/3"^3i/4^ + A2T4/3 1^3^vr2 zV^t:^) ] ^ ' 

is fixed uniquely by requiring the nonsingularity of the geometry at one order higher. In the 
above equation A denotes the overall energy scale defined by T = A • t^^/^ + . . ., tjq comes 
from the shear viscosity, while A^*^^ and are related to second order transport coefficients 
Ai and relaxation time rn (see [^]). The expansion in r~^/^ can be easily understood using 
hydrodynamics (see e.g. |^] for a detailed discussion), which is an effective description of 
systems at local equilibrium described by four degrees of freedom: the four-velocity 
{Ufj,u^ = —1) and the temperature T. Those quantities vary over space-time, however 
the scale set by their variations is much bigger than the microscopic scale given by the 
temperature. Thus the small parameter in this case is T~^V ^u'^ . In a boost-invariant setup 
the asymptotic scaling of the temperature obtained from the nonsingularity is t~^/^, the 
covariant derivative of the velocity^^ scales as which leads to the expansion parameter 



®The constant nn in first order, — — ''"I-" in the second order, . . . 

^"in a boost-invariant setup the fluid velocity takes the form — \ ■ [9r]'', but the metric in r 
coordinates has nonvanishing ChristofFel symbols leading to Vit — . 
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7--2/3 exact agreement with the one obtained on the gravity side. For phenomenological 
reasons (from causahty, see the discussion in e.g. ||2^) it is interesting to focus on an energy 
density expanded up to the second order in r~^/^ expansion, so that the geometry has to 
be solved up to the third order to fix ah the constants of interest. 

In |12| a surprising feature showed up: the pole contributions to the Kretschmann scalar 



at the third order cancelled, leaving a logarithmic singularity of the form log (3^/^ — v) . 
This singularity was present also in higher order curvature invariants in Feffer man- Graham 
variables, such as TZabcdT^'^^^^T^ef^^ ^ which led to the conjecture that the gravity dual 



to boost-invariant flow cannot be realized within the supergravity approximation 19 1. 
The resolution of this puzzle involves resummation of the series 

= nl {v) + -^nl {v) + -i^7^i (v) + -^7^i iv) + ... (s.u) 

which turns out to be equivalent to the change of variables from Feffer man- Graham to 



Eddington-Finkelstein metric [20|. Subsequently it was shown, that the geometry dual to 



the boost-invariant flow in the large r limit is well-defined at any order of r expansion 



1 21 ] and the usual black brane singularity is covered by an event horizon [gj, ^ . Thus the 
universal behavior of boost-invariant plasma at large proper-time (note that the energy 
density contains only one arbitrary constant, which sets the overall scale) possesses a 
regular bulk description, which could be constructed analytically in a perturbative manner. 
Moreover, as expected, it does not matter which coordinate frame in the bulk is used, so 
that in the rest of the paper all the calculations will be performed in Feffer man- Graham 
coordinates. 

4. Early time dynamics 
4.1 No scaling at early time 

The second regime where analytic methods can be applied is the domain of early times. 
At first sight it seems that the early-time problem could be addressed by the same method 
which was so useful for late time, namely solving the boost-invariant Einstein's equations 



( 4.10 ) without specifying a priori initial conditions. However after careful consideration 
a crucial difference turns out to occur: no scaling solution appears which would consist 
of a hierarchy of terms similar to ( p.5| ) and independent of the initial conditions. In fact, 
on a more physical ground, the initial conditions should play a crucial role at early times. 
This would indeed agree with the physical intuition which suggests that this regime should 
not be universal, but rather depend on initial conditions. However, as discussed in the 
following, the solutions are submitted nevertheless to nontrivial constraints if they are 
considered in the strong coupling limit. 

As will be shown in the rest of the paper the search for a scaling variable (as in Ref. 
||2^ ) is invalidated due to a subtlety. In order to see this one needs to recall the way how to 
identify the existence of a scaling variable. Consider the leading proper time dependence 
of the energy density in the regime of interest (r — > oo or r ^ 0) 

sir) ~ 1 (4.1) 
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and then solve the Einstein's equations ( |4.10D exactly as a power series around the boundary 
9,u = + g^J (r) + gl!J (r) / + gif) (^) / + . . . (4.2) 



The coefficients of this power series are exphcit s-dependent functions of proper time r 
and typicahy contain a couple of terms. If the part of gju) (t) which dominates at r ^ cxd 
is selected, the above power series becomes a power series of a scaling variable v = z/t^. 
If on the other hand the parts that dominate for r ^ are considered, it leads |23] to an 
expression of the form 

f(^^-) (4.3) 



for the metric coefficients. One consequently finds a unique solution in w for each s with a 
complex branch cut singularity for s > 0. This leads to the only possible value s = right 
at the margin of the allowed range, i.e. 

e(r) ~ const for r — > (4.4) 



out of the range of generic s |23|. 

However the situation when s = is special. Indeed, going back to the derivation of 
the early time scaling variable w = z/t, it turns out that the terms in g^^{T) which lead to 
it are all multiplied by a factor of s. Hence for s = 0, these terms vanish and a completely 
different hierarchy of terms appear. 

Indeed, if one performs for instance the power series expansion of a(r, z) starting from 
e(r) = l/r^ then the answer for the first three orders is 

~ z^ r-^ + z^ — -.-2 

\6 12 J 

+ z^ /_lr-2V--r-2^ + l/6T-2^s + -T-^-V-— r-^-vl + ... .(4.5) 
\ 16 6 ' 96 384 J 

Consider now the term proportional to z^ in (|4.5D . There are two structures z^/r^*"*"^-* and 
z^/r^*. The first of these leads to the early-time scaling variable proposed in |2^ since 
it dominates for nonzero s. However its coefficient is proportional to s so for s = it is 
absent and the only contribution comes from the second term which is not contained in 
the scaling-variable analysis. The same applies obviously already to the term in z^ /t^^'^'^^ 
as well as to all subsequent orders. 

Analysing the power series solutions in more detail, one finds that for a generic early 
time expansion of the energy density^^ 

00 

e(r) = 5]e2nr2", (4.6) 

n=0 

the coefficients a„ of the power series expansion in z of the metric coefficients at r = 

a(T = 0,z) = ^a„z4+2n ^4 7) 

n=0 



'^The restriction to even powers of r is discussed in one of the following sections. 
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depends on all coefficients e2n in ( f^.q) . Reformulating this observation, this means that 
each of the possible initial conditions ([4.7D leads to a distinct proper-time evolution (^), 
where the coefficients of the two power series are linked through the Einstein's equations. 
The following mapping is therefore obtained 



oo oo 

a(T = 0,z) = J^a„/+2- =^ e(r) =^e2nr2". (4.8) 

n=0 n=0 



On the other hand, certain nontrivial constraints limit the range of solutions. Among these, 
the nonsingularity condition on the metric-invariant properties of the geometry plays the 
crucial role. 

The derivation of this result together with the constraints including the nonsingularity 
argument restricting the early time expansion in terms of even powers of the proper time 
of the energy density, is the main result of this paper and is presented in section (|4.4| ) . 

The objective now is to analyze the qualitative properties of the evolution of the energy 
density starting from the initial conditions (fl.T]) . But before this can be done, the space 
of allowed initial conditions has to be investigated in more detail, as these have to satisfy 



a nonlinear constraint equation which is a part of the Einstein's equations ( 4.10 ). 
4.2 The FefFerman- Graham metric is singular at all times 

This section focuses on the restrictions on the initial conditions coming from Einstein's 
equations, which can be written in an equivalent simplified form as 

Rab + ^Gab = (4.9) 
or explicitly in Feffer man- Graham coordinates as 

(^rr) : b+2c-^{b+2c) + hp + 2c^)--{a-2b) = Sa"- — +(^--) {a'+b'+2c')\ , 
2 2 T \^ z \ 2 z J J 

{yy) ■.b-hb+\b-2a) + ]^{a+b+2c) {^+'^ = ("' + ^'+2^') 

(±±) ■.c-ac+^[a + b + 2c+'^ = e'^ |c" - ^ + (^^ - (a' + 6' + 2c')| , 
(rz) : 2b' + 4c' + b' ^6+ + 2cc' - a' ^6 + 2c + = 0, 

{zz) : a" + b" + 2c" - -{a' + b' + 2c') + -(a'^ + 6'^ + 2c'2) = 0. (4.10) 

z 2 

In the above expressions the ffi'st parenthesis is for the corresponding metric component, 
the dot for the proper-time r-derivative and the prime for the z-derivative. 

Expanding the {zz) component of Einstein's equations as a power series in z, and 
keeping only the leading terms for z ^ 0, one obtains the tracelessness condition — e (r) -|- 
P|l (r) -|- 2p_\_ (r) = of {T^u)- The {zz) equation thus encodes holographically the scale 
invariance of the boundary gauge theory. Moreover, that Einstein's equation does not 



- 12 - 



features any r-derivative. Hence, at each proper time, if the z-dependence of two of the 
functions a, b and c are given, the third one is completely determined by the (zz) equation 
and the requirement of a behavior near the boundary. 

However, for completeness, let us discuss also some of the other Einstein's equations. 
Expanding the (zt) component of Einstein's equations as a power series in z, one finds 
that the leading terms near the boundary give the relation + 2p_i + (e+p||)/T = which, 
assuming (T^) = 0, is equivalent to the energy conservation equation in the CFT. Hence, 
the {zt) and (zz) components of Einstein's equations together encodes holographically 
both energy conservation and conformal invariance in the CFT. 

One should remark that we have five equations ( 4.10| ) for only three unknown functions. 



There is indeed some redundancy in Einstein's equations ( [4.101 ). It is e.g. consistent to 
replace the first three equations (rr), (yy) and (-L-L) by the combination — (rr) + (yy) + 
2(_L_L) — e'^{zz), which writes 

c + 26 + 1^ = |26" + 4c" + b'^ + Sc'^ + 2b'c' - -\J - ^c'j . (4.11) 

Hence, the system of Einstein's equations is only of order one in r. The required initial 
conditions should thus be two of the three functions ao(z), 6o(-2) and co(z), where oq = 
oo {z) = a (r = 0, z) and analogous expressions for the others, because we have already at 
r = the constraint equation 

o!^ + 6^; + 2d^ + + 1(6^)2 + (c^)2 - i (a'o + + 2c'o) = 0, (4.12) 

We shall see in the next section that requiring a regular behavior of the solutions in the 
limit r — > further reduces the freedom for the initial condition to only one arbitrary 
function of z. 

The nonlinear character of equation ( 4.12| ) turns out to play a crucial role in the 



whole analysis. Naively, one would expect to be able to consider a small fluctuation over 
empty AdS'^ in the linearized approximation neglecting the quadratic terms in the above 
constraint. However it turns out, as will be shown below, that even if starting from e.g. 
an infinitesimal a, the nonlinear equation for c will always generate a singularity for some 
large but finite z\ Thus linearized fiuctuations cannot be used as initial conditions in the 
boost-invariant kinematics. Moreover, a singularity, presumably related to an horizon, has 
to be present from the outset already in the initial conditions. This point will be addressed 
later in the paper. 

To simplify (|4l2|) it is useful to introduce the following notation 



,(z2) = l6'o(z), 

w{z^) = ^c'o(z). (4.13) 

If the constraint equation has a regular solution, the derivatives Oq, 6q and Cq are bounded 
in the bulk and v^<^ and w vanish when z ^ 00. The constraint equation simplifies when 
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written in terms of ? and w 

+ + 2w' + + + 2?i;2 = 0, (4.14) 

where the prime now means the z^-derivative. Integrating the above equation assuming 
the regularity conditions from z = to oo gives 

/ {v' + + 2w')dz^ + I (t;2 + ^ 2w^)dz^ = 0. (4.15) 

The first integral vanishes due to the imposed boundary conditions leading to 

;>oo 

/ {v^ + ?^ + 2w'^)dz^ = 0. (4.16) 
Jo 

This equality however is satisfied only for v = q = w = thus the only regular solution 
is trivial - the vacuum AdS^. Therefore, at any time the metrics of interest must have 
a singularity at some value of z. In particular this will be the case even at r = 0, with 
the meaning that any nontrivial initial condition consistent with the Einstein's equations 
will lead to a metric singularity at some value of z. The nonsingularity constraint on the 
geometry requires that all the singularities apart from the one sitting at z = cx) will only 
be of coordinate nature. This provides a strong selection mechanism for the allowed initial 
conditions. 

4.3 Constraints at early time 

This section analyzes the constraints on initial conditions imposed by the Einstein's equa- 
tions together with the assumption of nonsingularity of the geometry at r = 0. The (rz) 
component equation (see ( [4.10| )) takes the form 

d,a-d,h = T ■{...), (4.17) 

which means that at r = 0, a(T, z) and 6(r, z) can differ only by a constant. Since both 
functions have to vanish at z = 0, this constant vanishes as well. Therefore 

a(0,z) = 6(0,z). (4.18) 



Incidentally the condition ( 4.18| ) ensures also that the geometry is nonsingular on the 
light-cone r = for generic z. 

In fact the constraint equation can be solved exactly leading to the space of solutions 
parametrized fully by the single function. The trick is to introduce the linear combinations 

v+ = —w — V, 

V- = w -V (4.19) 



for which the equation ( [4.14 ) becomes algebraic for v- (both and are understood as 



functions of z^). After trivial algebra one obtains 



2v'-vl. (4.20) 
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Therefore all solutions of the initial value nonlinear constraint equations are parametrized 
by an arbitrary function v-^-{z'^). The next step is to analyze what further conditions must 
be imposed on vji-{z^). First, since the the metric coefHcient functions have to vanish as z^ 
one gets 

~ ^£oz^ for z ~ 0. (4.21) 

Moreover, from the arguments in section 4.2 it follows that there is a singularity at some 
finite z = zq. The behavior at this singularity is constrained by demanding that this would 
be just a coordinate singularity and not a curvature blow-up. Assuming a power-like blow- 
up of v+ (z^) at z = zo, the regularity of the square of the Riemann tensor leads to the 
conclusion, that v+{z^) has to have a first order pole 

Zq z 

with residue 1. 

The coordinate singularity in at z = zq translates directly into the behavior of the 
metric coefficients around zq. This means that the proper time metric component has a 
second order zero at zq so that the metric at r = looks like 

ds^ = -l^(l-^\^ {dr^ + r^dy^} + ... + \dz^ (4.23) 

Z \ Zq J Z 

in the vicinity of z = zq, reminiscent of the behavior of a horizon in Fefferman-Graham 
coordinates. Note, however, that at r = the term in curly braces becomes dx^dx~ 
in contrast to a 'Schwarzschild' horizon where the corresponding structure is of the form 
— (1 — z/zQ)^dt^ + . . .. Presumably the metric singularity could be avoided by using a 
different set of coordinates (e.g. Eddington-Finkclstcin ones). Then it would be also 
possible to apply the framework of dynamical horizons to check and locate eventual trapped 
surfaces. However performing a change of variables from Fefferman-Graham to Eddington- 
Finkelstein ones is a formidable task, since contrary to late proper-time regimes, r ~ 
Fefferman-Graham is not mapped in the bulk into f ~ Eddington-Finkclstcin. This 
means that the change of variables cannot be performed perturbatively and is beyond the 
scope of analytical methods presented here. It would be however interesting to analyze 
in detail the local structure of the geometry sX z = zq including time-dependence. This 
problem is left for future work. 

It is interesting for further discussion to present explicit solutions of the constraints. 
One of them arises from the choice 

Vj^ = a(tana2:^ — tanhaz^). (4.24) 

The initial metric profiles may be integrated explicitly to obtain 

aQ{z) = 2 log cos az^, 

cq{z) = 2 log cosh az^. (4.25) 
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The above solution possesses a (coordinate - not leading to any curvature singularities) 
singularity at 

.0 = /J. (4.26) 
More generally one can parametrize in the following manner 

vA^') = l^oz!-^,f{z'), (4.27) 
where / (0) = 1, / [zq) = otherwise is a regular function of variable for z < oo. 



One has also to ascertain that v- obtained from (4.20), does not become complex 

Once the allowed initial conditions are under control one may proceed to find the 
solution of Einstein's equations starting with these initial data. 

4.4 Early time expansion of the energy density 

The Einstein's equations (|4.9|) can be solved for any energy density perturbatively in z^. 



Starting with some arbitrary energy density e, the first three nontrivial terms in the ex- 
pansion of a (r, z) warp factor take the form 



a (r, z) = -e(T) • z + 



At 12 j l6 ^ ^ 6 ^ ^ ^ ^ 16 ^ ^ 

128r^ 128r2 64r 384 ^ V ^ ^ 



This power series can be extended to an arbitrary order in by solving Einstein's 
equations and the only obstructions are of a purely computational nature. Generically 



terms in the expansion (4.28|) contain inverse powers of proper time multiplying the energy 



density and its derivatives. Assuming the energy density can be expanded in a regular 
power series around r = 0, the singular inverse powers of proper time in ( 4.28| ) will be 



present unless all the odd terms vanish. This requirement constrains the energy density in 
the early times domain to be a power series in even powers of r, namely 

e (r) = eo + ear^ + e^r"^ + ... . (4.29) 

On the other hand, taking r to be zero in equation ( [4.28| ) gives the relation between the 
early times energy density and the profile of qq (z), which has been signaled in subsection 
( |4.1| ). Basically expanding the initial profile of the metric in the radial AdS variable near 
the boundary and comparing with (|]2|) at r = 

a (0, z) = ao (z) = -e^ ■ z^ - '^€2 ■ z^ + (-'-^ - f \ ■ z' + . . . (4.30) 



allows one to solve for all the e2i sitting in ( 4.29 ). This pattern continues to any order 



in z (and thus r) expansion. Note also, that each parameter in the expansion of energy 
density around r = is an independent dimensionful quantity and in order to fully specify 
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the initial conditions infinitely many such terms are needed. This is in stark contrast to 
the late time behavior, where only one dimensionful constant (A) appears. Physically, this 
is in agreement with the fact that in the thermally equilibrated final stages of Bjorken 
expansion all differences due to initial data should have been washed away by dissipative 
effects and the only parameter characterizing the flow is an overall energy scale (given e.g. 
as the energy density at a certain fixed proper time tq). 

The strategy for finding the time evolution from given initial data is to use the Ein- 
stein's equations to generate the expansion of the metric warp-factors for sufficiently high 
order (the bigger, the better, but for technical reasons this has been achieved up to the 
order or equivalently t^^ in most cases) and then for given regular (with finite curvature 
evaluated on initial data for z < oo) initial profile generate the power series for the energy 
density at early times. Thus the AdS/CFT correspondence beautifully sets the allowed 
initial conditions for gauge theory dynamics using the gravitational description. However 
since the series ( [4.29 ) has generically a finite radius of convergence, some resummation 
method is needed in order to extend it to larger proper times. 



5. Transition to the hydrodynamic regime 
5.1 Resummation scheme for the energy density 

It is desirable now to find a suitable numerical approach, relating the solutions of the 
Einstein's equations ( [4. 101 ) to the physical quantities, energy density and pressures, as a 
function of r. 

Though early time dynamics of the field theory is dictated by the initial conditions, 
after certain time the system is expected to settle down to local equilibrium. For r suf- 
ficiently large the plasma should exhibit the universal hydrodynamic behavior, where the 
only trace of the initial conditions is given by the overall scale A, see ( 3.13| ). In order to 



track the dynamics of the system with sufficient accuracy from r = to r ^> 1, numerical 
methods are needed. In particular the early-times power series for the energy density has 
a finite radius of convergence and a resummation is needed in order to find its behavior for 
large r. 

We would like to address the following questions. Firstly, whether during the evolution 
from some (generic) initial data at r = one can observe a passage to the asymptotic perfect 
fluid behaviour 

6(r)~-i^ + ... (5.1) 
To be more general one may try to determine the asymptotic exponent s in 

e(T)~l + ..., (5.2) 

and determine whether it is significantly different (as is expected) from the free streaming 
value s = 1. 
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Figure 2: Approximate value of s obtained from the logarithmic derivative and Pade approximation 
for A) Ncut = 32 (dotted line), B) Ncut = 40 (dashed line) and C) Ncut = 48 (solid line) for initial 
profile v^"^ i^'^)- Two horizontal lines denote s = 1 (free streaming scenario) and s = 4/3 (perfect 
fluid case). 



As explained above, it is quite difficult to answer this question as the energy density 
in the early time regime takes the form 

6 (r) = 60 + esr^ + . . . + esTv.^.r^^™* + . . . (5.3) 

where Ncut is a natural number denoting the cut-off up to which the evaluation of the 
energy density from the initial profile in the bulk has been performed. In most cases 
Ncut = 40 and increasing this accuracy is difficult. Moreover the series has a finite range 



of convergence. In order to estimate the asymptotic exponent s appearing in (5.2) it is 
convenient to express it through a logarithmic derivative 

s = -lim^^ooT • ^loge(T). (5.4) 



and perform a Pade approximation (of order {Ncut, Ncut)) to the r.h.s. of (5.4) with e(r) 
substituted with our power series ( |5.3D : 

(2)^2 _L _L A^^'=nt)^2N^ut 

(5.5) 



s = 4/3 corresponds then to the perfect fluid case and s = 1 to the free-streaming scenario. 
Of course it is not expected that ( |5.5D will give s = 4/3 exactly^^. 

Let us consider the three following examples of the initial profile, satisfying the curva- 
ture nonsingularity constraint, the first one being ( |4.24| ), the second its deformed variant. 



and the third one from the family (4.27): 



(1) 

+ 




= tan 








= tan 




(3) 




2 




+ 




3 1 


— z 



, / 2 ^ 

tanh { z -\ 

V 6 

<M-^) = ^T— 2{l + i-^V (5-6) 



^^Such a Pade approximation has, by construction, a different subleading large r behaviour from viscous 
hydrodynamics . 
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The obtained values of s for the first profile, for which we had a power series in r 
up to order r^°°, are Sapprox = 1.1667, 1.1923, 1.2025, 1.2025 and 1.2101 respectively for 
^cut = 32, 36, 40, 44 and 48. The approximate value of s is closer to s = 4/3 than to 
s = 1 for largest Ncut- For the second and third profile the results are inconclusive, since 
the corresponding energy densities are provided with worse accuracy. 

In the rest of the text, we shall keep the asymptotic perfect fluid value s = 4/3. Then, 
a suitable Fade resummation scheme, with the required asymptotia can be given by 

M (2) 2 , , (Ncut~2) N,^^^2 

— 4°) + egv^ + . . . + 4^--2),^...+2 • 

where both ej^ and ef} are obtained by expanding (^) around r = and comparing 
with ( [5.31) . Such a resummation imposes the correct asymptotic behavior, but differs, by 
construction, in the subleading behavior with viscous hydrodynamics. 

^ + (5.8) 



^4/3 ^2 

(the first subleading piece scales as whereas the correct scaling is r^^/'^). This difference 
is not substantial and can be cured by more refined resummation schemes^^. 

Despite its simplicity, it turns out that the resummation ( ^.7] ) works pretty well (the 
results seems to converge well with increasing cut-off, see Figs. 2 and 4 in extending the 
energy density beyond the convergence radius of the early-time power series (which is its 
main task) and in providing a qualitative picture of dynamics. 

5.2 Qualitative features of the approach to local equilibrium 



The method (|5.7| ) can be used to study the approximate behavior of the energy density as 
a function of time for large enough times to see local equilibration. Fig. 3 shows the plots 
of energy density as a function of proper time for the three profiles (5.6) obtained for the 



highest cut-offs. The results seems to converge well for the first profile, see Fig. 4. Energy 
densities obtained for these profiles differ at the initial stages, whereas in the late-time 
regimes both seem to approach local equilibrium. A measure of the local equilibrium is the 
relative difference between the transverse and perpendicular pressures defined as 

Ap(r) = 1-^44 (5.9) 
P± (t) 

When this quantity is close to zero, it signals the isotropisation indicating local equi- 
librium, while a value of order one is an indication in favor of the free streaming scenario. 
Fig 5 shows the plot of the relative difference of pressures ( |5.9D . It is interesting to note 
the rapid fall-off of the pressure difference on a scale r = 0(1). This fall-off appears to 
be stable after different numerical checks. Interestingly enough, there is a bump which 
prevents the pressure to reach isotropy before r = 0(5). However, the Fade approximants 



^^An other possible issue are roots of the denominator lying within the range (0, oo). If such feature is 
encountered, it can be interpreted as an artifact of approximation without real physical significance. 
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A) B) C) 

Figure 3: A) Energy density ei (r) as a function of proper-time r obtained from Pade approxi- 
mation for cut-off Ncut — 46 and initial profile w^^"* (z^) in the bulk; B) Energy density £2 (t) for 
the second profile. C) Energy density £3 (r) as a function of proper-time r obtained from Pade 
approximation for cut-off N^ut = 34 and initial profile vf in the bulk. 




Figure 4: Relative difference between the energy densities for the first profile for cut-offs Ncut = 16 
and Ncut — 46 does not exceed 10%. 



for the pressure difference are less stable than for the energy density and the differences 
appear after the bump. The second profile which is a slight deformation of the first one 
does not seem to exhibit this bump. In any case, it would be physically interesting to 
check whether this phenomenon of rapid fall-off but incomplete isotropization is or not 
a characteristic feature of the strong coupling evolution. We intend to analyze this issue 
using numerical methods for solving Einstein's equations in future work. 

6. Summary 

This paper studies the early-times dynamics of boost-invariant plasma using analytical 
methods. The motivations for this work are both phenomenological and theoretical. On the 
one hand, the recent findings of the RHIC experiment suggested that locally equilibrated 
nuclear matter behaves as an almost perfect fluid, which presumably indicates a strongly 
coupled regime (at least at this stage of its evolution) of the underlying gauge field theory, 
i.e. QCD. Despite the fact that realistic dynamics of the QCD after the collision would 
require the understanding of rapidity dependence [i.e. deviations from boost-invariant dy- 
namics, which is only approximately valid in the central rapidity region) and might be 
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Figure 5: Relative difference in pressures for the first (left) and second (right) profiles - for r « f 
there is a rapid fall-off but (perhaps, see text) does not reach yet a complete izotropization for the 
first profile. 

driven as well by the perturbative effects (as for the Color Glass Condensate initial 
conditions), it is still interesting to consider the boost-invariant expansion of a conformal 
plasma at strong coupling as a useful toy-model for e.g. estimates of the thermalization 
time. On the other hand, theoretically, the boost-invariant plasma evolution is a workable 
example of a dynamical system at strong coupling, which is interesting on its own. The 
modern developments translates (using the AdS/CFT correspondence) the dynamics of 
the strongly coupled gauge theories into the evolution of higher dimensional space-time 
equipped with a nontrivial metric. Thus, as suggested previously by various authors (e.g. 
pO|]) the thermalization of the excited gauge theory matter should be dual to black hole 
(or black brane) formation (see ^] for concrete realizations of this observation) , which 
is obtained in late time as the dual of the Bjorken hydrodynamical flow. This subject is 
very fresh and thus any work which may shed light on this fascinating process is valuable. 

The main result of the studies presented here is that the boost-invariant dynamics of 
a strongly coupled conformal plasma is sensitive to the initial conditions. This contradicts 
the scaling hypothesis |2^, which, analogous to the late-time case Q would indicate some 
uniqueness of the early-times solution. In fact the scaling does not occur due to a subtlety 
- the a priori dominant scaling contributions vanish precisely in the limit s — > 0. (see 
subsection ( [4.1D for details). The correct physical picture leads to a link between the 
early-time expansion of the energy density and the initial profile of the bulk metric. We 
find a quite general result that at all times, including the initial one, a singularity of at 
least one metric coefficient should develop in the bulk. Hence, the requirement that this 
does not lead to real curvature singularities already at initial time is a basic constraint 
on the possible initial conditions. The analysis of the possible curvature singularities in 
the initial data at finite distance from the boundary fixed the early times power series 
for the energy density to contain only even powers of proper time. We have shown that 
solving the nonlinear constraint equation in the Fefferman - Graham coordinates leads to 
the conclusion that the initial data must contain a coordinate singularity in the bulk of 
AdS which may signal the presence of a dynamical horizon (we are making a distinction 
w.r.t. the event horizon which is a global notion) right from the start of the evolution. 
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This suggests that the process of thermahzation in this boost-invariant setting does not 
amount to a formation of a black hole, as was suggested up till now, but rather amounts 
to a qualitative difference in the behaviour of the horizon. We plan to investigate these 
issues in more detail in future work using numerical methods. 

Interesting further directions of study include numerical investigation of the bulk evo- 
lution dual to the boost-invariant flow. This can be achieved using the methods presented 
in 1 31] and may provide precise numerical results covering both early, intermediate and late 
time regimes (see |33| for an interesting study of boost invariant flow sourced by bound- 
ary metric perturbations). It is important to understand more qualitatively the relation 
between the initial conditions in the bulk and the shape of the energy density as function 
of proper time, which would give the more precise estimates on the thermahzation time. 
Moreover the methods developed in this paper are well suited to reconsider the problem of 
plasma isotropisation posed in [^] using analytical methods (note that anisotropic energy- 
momentum tensor should reach the equilibrium exponentially fast, whereas hydrodynamic 
evolution leaves a power- like tail). Finally the most interesting, yet highly non-trivial, ex- 
tensions of the AdS / CFT program for the dynamical evolution of a plasma are the studies 
of the initial conditions such as shock waves collision using the dual gravity picture (see for 
instance 35 1 for some preliminary attempts). We hope that the general properties and 
their numerical implementation we found can serve as a testing ground for the investigation 
on the proper initial conditions and evolution of the plasma. The authors plan to address 
some of the mentioned issues in the nearest future. 
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